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future ICASE reports in this area of research will be printed with a green cover. Applied 
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THE NONLINEAR EVOLUTION OF MODES ON UNSTABLE 

STRATIFIED SHEAR LAYERS 


Nicholas Blackaby, Andrew Dando & Philip Halt 1 
Department of Mathematics 
Oxford Road 
University of Manchester 
Manchester, M13 9PL 
UNITED KINGDOM 

ABSTRACT 

The nonlinear development of disturbances in stratified shear flows (having a local Richard- 
son number of value less than one quarter) is considered. Such modes are initially fast grow- 
ing but, like related studies, we assume that the viscous, non-parallel spreading of the shear 
layer results in them evolving in a linear fashion until they reach a position where their 
amplitudes are large enough and their growth rates have diminished sufficiently so that am- 
plitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer 
theories. Four different basic integro-differential amplitude equations are possible, including 
one due to a novel mechanism; the relevant choice of amplitude equation, at a particular 
instance, being dependent on the relative sizes of the disturbance amplitude, the growth rate 
of the disturbance, its wavenumber and the viscosity of the fluid. This richness of choice of 
possible nonlinearities arises mathematically from the indicial Frobenius roots of the govern- 
ing linear inviscid equation (the Taylor-Goldstein equation) not, in general, differing by an 
integer. The initial nonlinear evolution of a mode will be governed by an integro-differential 
amplitude equations with a cubic, nonlinearity but the resulting significant increase in the 
size of the disturbance’s amplitude leads on to the next stage of the evolution process where 
the evolution of the mode is governed by an integro-differential amplitude equations with a 
quintic nonlinearity. Continued growth of the disturbance amplitude is expected during this 
stage, resulting in the effects of nonlinearity spreading to outside the critical level, by which 
time the flow has become fully nonlinear. 


'This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-19480 while the three authors were in residence at the Institute for Computer Applications 
in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23681. 
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1 Introduction 


This paper extends the theories of Churilov & Shukhman (1988) and Goldstein & Leib 
(1989) to the case of (non-marginally) unstable stratified flows. Such flows, in which the 
vertical density variation is insufficient to overcome the strongly destabilising effects of the 
vertical velocity variation, occur naturally in many situations. The mathematical formulation 
of this problem has great similarity with, amongst others, that for the related problem of 
buoyancy-driven instabilities above a heated plate and with the, at first sight, less physically 
related problem of vortex disturbances driven by wall curvature occurring in 3D boundary 
layers. Thus, whilst here we concentrate solely on a particular model stratified shear flow, the 
mathematical theory has much wider applications to both geophysical and aerodynamical flow 
situations. 

The stratification of a shear flow is usually characterised by a physical parameter called the 
Richardson number (herein after denoted J); it is a measure of the ratio of vertical density 
variation to the scale of vertical velocity shear. The Richardson number can be both positive 
and negative; the former case corresponding to lighter fluid lying above heavier fluid, and the 
latter case corresponding to the lighter fluid lying below. The stability of a stratified shear flow 
was first considered independently by Taylor (1931) and Goldstein (1931). They concluded that 
a multi-layer system of homogeneous fluids can not be used to approximate the stability of a 
heterogeneous fluid, and that for the flow models considered, the flow is stable for Richardson 
numbers greater than one quarter. 

The flow model adopted in this paper was introduced by Drazin (1958); the vertical ve- 
locity and density distribution are represented by hyperbolic tangent and exponential profiles, 
respectively. This model enables an exact solution of the governing linear stability equation, 
the Taylor-Goldstein equation, for neutral disturbances and the same upper bound on the 
Richardson number is found. Soon after Drazin’s paper appeared, Miles (1961) presented 
several theorems concerning properties of solutions to the Taylor-Goldstein equation. In par- 
ticular, he proved the two important results that in general flows are stable for J > 1/4 and 
that neutral modes are proportional to just one of the associated Frobenius series near the 
critical level. The continuation of the linear eigenfunctions in the neighbourhood of the critical 
levels was considered by Koppel (1964). 

Several linear results for the alternative Holmboe model were computed by Hazel (1972), 
for a variety of boundary conditions; however, our interest herein is in the nonlinear critical- 
layer theory and not with a discussion about which model is superior. Around the same time, 
nonlinear studies started to appear eg. Maslowe (1972) using the equilibrium critical layer 
theory of Benney & Bergeron (1969). Subsequent nonlinear studies were also mainly for the 
marginal instability case (0 < 1/4 — J 1) eg. by Maslowe (1977) who considered a viscous 
critical layer, and by Brown & Stewartson (1978) who considered a time-dependent critical 
layer for small time. These two theories were unified by Brown, Rosen &; Maslowe (1981). The 
latter paper was corrected to some extent by Churilov & Shukhman (1987); soon after, Churilov 
& Shukhman (1988) essentially extended and corrected the work of Brown & Stewartson (1978) 
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and Churilov k Shukhman (1987). Recently, Troitskaya (1991) has considered nonlinear effects 
for J > 1/4. 

There have been many more studies of (non-stratified) shear layers, in which there is no 
vertical density variation. The linear stability of such flows is, in general, governed by the 
familiar Rayleigh equation, to which the Taylor-Goldstein equation reduces for zero Richard- 
son number. The nonlinear stability of such flows has received much attention over the past 
couple of decades. Benney k Bergeron (1969) developed the so-called equilibrium critical layer 
theory; here the mode is treated as ‘quasi-steady’ inside the critical layer as well as outside 
it. Nonlinearity affects the jump imposed across the critical layer and hence leads to modified 
results for the neutral (equilibrium) modes. Haberman (1972) extended the theory to include 
critical layers where viscosity is also significant. The first studies of non-equilibrium critical 
layers include the papers by Stewartson (1978), Warn k Warn (1978), Hickernell (1984) and 
Gajjar k Smith (1985). The key paper by Hickernell (1984) concerned a shear layer affected 
by Coriolis (rotational) effects; here the weakly nonlinear theory leads to an integro-differential 
equation rather th^ti* the (previously) more familiar Landau equation with ‘polynomial’ non- 
linear terms. In fact, such integro-differential equations result naturally from non-equilibrium 
nonlinear critical layer theories when the shear layer is coupled with other physical factors such 
as, for instance, Coriolis effects (eg. Hickernell, 1984; Shukhman, 1991); compressibility effects 
(eg. Goldstein k Leib, 1989); three-dimensionality effects (eg. Goldstein k Choi, 1989; Wu 
et al, 1993); and buoyancy effects (eg. Churilov k Shukhman, 1988). Moreover, in a related 
paper (Blackaby, Dando k Hall, 1993), we demonstrate that such an integro-differential type 
of equation can be derived to describe the nonlinear evolution of the inviscid Gortler modes 
studied by Bassom k Hall (1991). However, the case of a ‘simple’ shear layer, not affected 
by any additional physical factors, is a special case in the sense that it does not lead to an 
integro-differential equation; instead, Goldstein k Leib (1988), found that the nonlinear evolu- 
tion of a disturbance was governed by the full unsteady nonlinear critical-layer equations. This 
difference is due to the additional physical factors, of the former cases, resulting in stronger 
singularities of the inviscid disturbance quantities at the critical level. 

At first sight, it appears that weakly nonlinear theories can only be usefully applied to 
marginally unstable flows; they rely on small growth rates and so the unstable disturbance of 
concern must be near to a neutral state. Thus, it was believed that such theories are of no 
use in describing the initial evolution of ‘far-from-neutral’ unstable modes. However, several 
recent studies have derived integro-differential equations, using weakly nonlinear theories, to 
describe the nonlinear evolution of (general) unstable modes on a variety of shear layers (see the 
previous paragraph). These studies are based on the assumption that, in actual physical flow 
situations, shear layer spreading or other external changes/effects would result in the otherwise 
relatively unstable modes having their growth rates diminished in real terms, so that a weakly- 
nonlinear critical-layer theory becomes appropriate; in fact, we adopt the same argument for 
this present study. The general proposal/argument is supported by the findings of Michalke 
(1964), Crighton k Gaster (1976) and Hultgren (1992). 

In this study, we use the unsteady critical-layer and weakly-nonlinear theories to describe 
the nonlinear development and evolution of unstable linear disturbance modes on stratified 
shear flows where the Richardson number J takes values less than one quarter. In particular, 
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we consider the temporal evolution of two-dimensional (2D) modes on Drazm’s (1958) model 
flow; we note however that the analysis can easily be modified for the spatial evolution case 
and/or a different choice of model flow. In addition to not complicating the analysis, the choice 
of considering only 2D disturbances is justified by noting (i) the model flow is also only 2D, and 
(ii), that Squire’s theorem holds for such flows where Boussinesq’s approximation is used (see 
Koppel, 1964). Perhaps the main assumption made, in this first such study for J < 1/4, is that 
the Reynolds number of the flow is large enough such that the critical layer is not significantly 
affected by viscous effects to the orders considered. Such an assumption appears justified based 
on our results. 

Although this paper is solely concerned with the nonlinear evolution of (non-marginally) 
unstable stratified shear flow which has important geophysical applications, the actual initial 
motivation for the study was the authors’ desire to develop a theory to describe the nonlinear 
evolution of the inviscid Gortler modes considered initially by Bassom & Hall (1991), and lately 
by Dando (1992), Blackaby k Choudhari (1993). The first two of the above papers demonstrate 
that in the presence of a relatively weak cross flow, longitudinal vortex disturbances of all 
wavelengths are stabilised such that the inviscid modes possess some of the largest growth rates 
whilst also being neutral at' certain other wavenumbers. Their governing equation is similar 
to the Taylor-Goldstein equation which, as mentioned earlier, governs the linear stability of 
stratified shear flows. In fact Blackaby k Choudhari (1993) have illustrated the close connection 
between the two problems and propose a definition of a generalised Richardson number for 
such centrifically-driven instabilities. The ideas developed here in this paper also have obvious 
applications to inviscid modes in flow above a heated plate, similar to those considered by Hall 
k Morris (1992). 

Whilst different from the approach adopted in this study, there are alternate/ complementary 
nonlinear theories that have been developed recently in which two or more of the flow distur- 
bances mutually interact. Such theories generally require smaller disturbance amplitudes but 
may also need the disturbances to exist in specific ‘configurations’. These other theories are 
generally referred to as wave/ wave and vortex/ wave interactions. For a discussion of wave/ wave 
interactions and resonant-triads the reader is directed to the book by Craik (1985) and the pa- 
pers by Goldstein k Lee (1992), Wu (1992). Strongly nonlinear vortex/wave interactions were 
first looked at by Hall k Smith (1991) and their ideas were clarified and extended by Brown, 
Brown, Smith k Timoshin (1993). Dando (1993) has looked at this type of interaction in a 
heated boundary layer where both streamwise vortices and inviscid travelling waves are present. 

The format of the rest of the paper is as follows. In the next section we present some 
background details of the flow concerning us in this paper, namely the stratified shear flow 
model due to Drazin (1958), as well as introducing some notation and concepts ready for the 
proceeding sections. In §3 the flow outside the critical layer is considered, whilst in §4 we 
consider the flow inside the critical layer. In §5 we discuss the range of validity and solution 
properties of the evolution equations. Finally, in §6 we draw some conclusions. 
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2 The stratified shear flow model 


Following Churilov & Shukhman (1987,88) we choose to consider the simple but realistic 
stratified shear flow model of Drazin (1958), where the unperturbed velocity in the x direction, 
uq, and the density, po, are dependent on the vertical coordinate y (in fact u 0 = tanh y and 
Po = Poc — y). The governing equations of motion in the Boussinesq approximation can then 
be written in the form 

- J-^ + {AV>,V>} = 

+ {p, 0} = (2.1a, 6) 

where {a, b} = [da/ dx){dbj dy) — ( da/dy)(db/dx ); 0 is the stream function; p is the density; rj 
the Prandtl number; k the thermometric conductivity and J the Richardson number. 

The shape of the stream function of a two-dimensional disturbance mode of the inviscid 
linear problem, denoted ipi say, is found to satisfy a Taylor-Goldstein equation of the form 

UtP = 0, (2.2) 


where the operators L; have the form 

T _ d 2 ( 2 _ 2 sinh ysech 3 y J \ 

1 ~ dy 2 \ (tanhy — c) (tanh y — c) 2 J 

Here k and c are respectively the wavenumber and wavespeed of the infinitesimal distur- 
bance, whilst the index / corresponds to the harmonic being considered (see later). Note that 
the L; are singular where tanh y — c\ at such locations, usually referred to as critical levels, the 
assumptions employed in the linear inviscid approach are no longer valid and special attention 
must be paid. Thin regions around such levels, usually referred to as critical layers, must be 
introduced into the mathematical model; such layers are generally the first to ‘feel’ the effects 
of increasing disturbance amplitude. 

Drazin (1958) found neutral eigensolutions of the form 

= £±*My), 0a(y) = sinh 1-fc2 |y| sech y , 

with c = 0 and J = A: 2 (l — fc 2 ), (2.4a — d ) 

satisfying (2.2) coupled with the boundary conditions — » 0 as y — ► ±oo. Here the ± on B 
are related to vertical position with respect to the critical level at y = 0 i.e. y > 0 and y < 0 
respectively. The relationship between the amplitudes B + and £_ will follow from matching 
to the inner problem (cf. Miles, 1961). In fact, in §4 we shall find that J9_ = i~ l ~ 2l/ B where 
we have written B = B + and the Frobenius root v is defined below. 

For later convenience, we introduce the quantity 

!/ = (2.5) 
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note that, unlike Miles (1961), we choose to include the “±” in our actual definition of v — 
this is also done for convenience. It easily follows from (2.4,5) that for Drazin’s model flow 
u = 1/2 - k 2 and thus the appropriate root is immediately determined. However, for general 
flows, an analytical solution is not possible and the relevant choice of sign for given k, c and J 
can only be determined by inspecting the shape of the eigenfunction obtained from a numerical 
solution; however, there does appear to be a common pattern between location on neutral curve 
and which sign in (2.5) is appropriate (see Blackaby k Choudhari, 1993). Note that half of 
Drazin’s neutral curve corresponds to the positive choice of sign and the other ‘half to the 
negative choice (see Figure 2.1); moreover, note that Drazin s solution is also valid for J < 0 
which corresponds physically to buoyancy driven flows. For J < 0, it appears that the — 
option is always relevant; thus the eigenfunction also is singular at the critical level. 

A study of the fundamental and other low harmonics outside and inside the critical layer 
is necessary to derive the desired amplitude equations; in the next section we consider the 
flow outside of the critical layer. Here the details are dependent on the flow model under 
consideration, but the method is quite general and can be applied to other flows. In §4 we shall 
see that the critical-layer analysis is almost entirely independent of the stratified shear flow 
being considered. 


3 Outside the critical layer 

We consider a wave of sufficiently small amplitude in the neighbourhood of a general point 
(J,k) on the neutral curve of Drazin (1954) (see Figure 2.1). Note that, since we are not 
considering the initial evolution of a mode on a marginally unstable flow, we are not restricted 
to the sole case J ~ 1/4. In fact, we begin by allowing J to take any value < 1/4 (i.e. we choose 
to exclude the special case J = 1/4), and then see which values of v = v(J,k) require special 
attention later. Note that J ^ 1/4 for all but one point on each of Drazin’s and Hazels’ (1972) 
neutral curves; moreover, the generalised Richardson number is always negative for neutral 
inviscid longitudinal vortices in 3D boundary layers. In fact, the theory of this and the next 
section can be regarded as the extension of the work of Churilov k Shukhman (1988) to the 
case v ^ 0, i.e. to a far wider range of problems. At the outset, we also choose to exclude the 
special case J — 0, which corresponds to a flow with no stratification (this case is considered 
by Goldstein k Leib, 1988). 

Following Churilov k Shukhman (1988), we introduce the small parameter e, characterising 
the magnitude of the mode. It is also necessary to introduce the ‘slow’ evolutionary time t = (it, 
with (i also small i.e. the amplitudes B± appearing in (2.4a) are considered to be functions of 
r. Note that although the time-scale r is slow with respect to the inviscid timescale t, it is 
still much ‘faster’ than the timescale t ~ /c -1 of the viscous spreading of the base shear flow, 
provided 

k < (i (< 1 ); 

this restriction is not severe due to the large size of the Reynolds number Re ~ k~ 1 for flows of 
practical interest. Later we shall have to relate the small parameters k, e and (i to one another 
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so that we can balance the linear terms of the amplitude equation with possible nonlinear ones. 
We write 

J — Jo + pJ\, (3-1) 

where Jo corresponds to a general point on the neutral curve and J\ < 0 is order one. The 
streamfunction and density are expanded as Fourier series in x, 

00 00 

v> = ln(cosh y) + ^ 4>,(t, y)e ,klx , p = p 0c -y + M r >y) e ‘ Wx > (3.2a, 6) 


l=—oo 


l —— c 


with boundary conditions 

ipi, Pi 0 as 3/ — ► ±oo. (3.3) 

In the present study, we need to consider the fundamental, the zeroth harmonic, the second 
harmonic and the third harmonic. Equations (2.1a,b) are solved both outside and inside the 
critical layer and then these solutions are matched to obtain the evolution equations. Outside 
the critical layer the fundamental harmonic dominates the perturbation and has an amplitude 
of order e. The other harmonics are the result of self-interactions, with the zeroth and second 
having amplitudes of order e 2 and the third of order e 3 . 


3.1 The fundamental harmonic 

In order to derive the evolution equation, the following terms of the expansion 

V>i = + e/i^ 2) + ck^ 3) + • • • , (3.4) 

need to be considered. Here ijj{^ is the neutral mode of the inviscid linear problem; takes 
into account the T-dependence of the solution; and is a correction to for dissipative 
(viscous) effects. The corresponding pi expansion includes analogous terms. Thus 


V’l 1 * = B±(T)i>a(y) where rp a (y) = sinh* +l ' |y|sechy; 


i£+* 


we note that 

V>i a) ~ B ± (T)\y\*^ , as y 
The second term in the expansion (3.4) satisfies 

LiV>i 2) = Qu i’? -» 0 as y 
J\B±ip a 2t0 a dB± 


0 . 


where 


Q i = 


±oo, 

Jo 


\ cosh 2 y tanh 2 y J 


(3.5 a, 6) 
(3.6) 

(3.7a) 

(3.76) 


tanh 2 y k tanh y dr \cosh 2 j/ tanh 2 j/ / 

The solution can be considered to be the sum, ip[ 2 ^ — V’ip/ + V’icf* a particular integral 
of (3.7), ipipj say, and the complementary function, ip[cF sa y- As y — ► 0, it follows from Taylor 
(1931) that 

-&V = B ± a$\y\i- l -(l + 0(|j|-')) + B ± b^\y\l-{l + Odsfi" 1 )), (3.8) 
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where and 6 jj! are constants as yet undetermined. Note that if the Frobenius roots j i [t'l 
differ by an integer then (3.8) is no longer appropriate (logarithms are needed). As such cases 
(i/ = |m; m integer) are isolated, we choose not to concern ourselves with them (and their 
immediate neighbourhood) in this paper. 

A solvability condition for the above boundary-value problem (3.7) is required. Note that: 
(i) the operator Lj is self-adjoint away from the critical level y = 0, and (ii) the right-hand side 
is singular at y = 0. Rather than closely follow the method of Churilov & Shukhman (1987) 
for deriving the modified solvability condition, for Jo 7^ 1/4 R i s more convenient to adopt the 
related approach employed by, for instance, Hickernell (1984). 

The solvability condition is derived by multiplying both sides of equation (3.7a) by V’o and 
integrating over all y, excluding the (sole) critical layer at y = 0. After integrating by^parts; 
imposing the boundary conditions at y = ±00; and the asymptotic forms of ift a and fp\cF 35 
y — *■ 0, we find that 


-j V’aQi dy = - [t/’icWl “ V>i cWa] 0 _ = -~g { B l b ?+ + > ( 3 - 9a ) 

where the bar through the integral sign indicates that the finite part of the integral should be 
taken. After substituting for Qi and using the relations B- — i~ l ~ 2v B, B + = B (to be derived 
in the next section), the solvability condition becomes 


4i 1+2l/ Ji 


cos(7tj/) 


dB Jihk . 


= -2vB{b { ?l - r^b?}), 


(3.96) 


where 


h = { 

Jo 


0O sinh 2l/ y 


Jo 


fOO 


coshy ^cosh 2 y tanh 2 y I * Jo 


dy and I 2 = 4 sinh 2,/ 1 ydy. 


(3.9c, d) 


Four relations involving 6$ and b{ 2 } are determined from the inner problem considered in 
the next section, thus determining the possible evolution equations for the wave amplitude 
B(r). 


3.2 The zeroth, second and third harmonics 

The presence of these terms is due to the process of harmonic generation, i.e. due to 
nonlinearity the fundamental of 0(c) generates the zeroth and second of 0(c 2 ) and so on. The 
zeroth and second harmonics are expanded in the form 


*, = eVS” + £V4 2) + <^4 3) + • • ■ , (3-10) 

<h = ‘VS 11 + *V4 2) + = 2 *-/4 3) + • • • , (3-ii) 


and sim 
namely 


ilarly for p 0 and p 2 . It is only necessary to consider the first term of the third harmonic, 


V>3 = + 


(3.12) 
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and the similarly for p 3 . The resulting sets of equations from expansions (3.10-12) do not have 
simple exact solutions when Jo ^ 1/4; however it is, in general, a relatively simple process, 
if a little tedious, to deduce their asymptotic forms as the critical level is approached. The 
significant terms of these asymptotes are quoted in the next subsection. 

However, the first term of the second harmonic requires extra attention as it can lead to the 
largest nonlinearity in certain circumstances when v is positive (see §4). The necessary analysis 
is simply the generalisation of that for the case v = 0 (J = 1/4) presented by Churilov h 
Shukhman (1987) in their unpublished Appendix B. Since this important work is unfortunately 
not readily available for reference, we take the liberty of presenting, in Appendix A of this 
paper, more details of their method than we would otherwise. 

3.3 The asymptotic expansions as y — ► 0 

In terms of the new variable Y = p~ l y and the functions 

^ ^ - (1 + (3.13) 

where pi = —2 dpi/dy, the asymptotes for the fundamental, second, zeroth and third harmonics 
as y — ► 0 are 




B±\Y\*+ V - 


*+„ i(l + 2u) dB ± \Y\* +v 


+e/i» +t/ 


2k dr Y 


a [ ?lB ± - 4^ In 


+enp 


,± ~ I 21/ 


i/^i] in 


l!+" 




+ u 


48k 

+d ± |y|^ + ..., 


{ v (5-2u) + (l-2v)}B ± !-jz r + 


<f> i = cpi* 1 '. 0 + ep* - *' 


-4vB ± b%l 


(1-2 u) 






Z^lE±-\y\k+ 


i/(l -j- 2 1/) 


+£K^”2 + * / 


i(T) - 1)(1 - 4i/ a )B± \Y\ 


8k 


Y 3 


+ 


(3.14a, b) 


- .2 (1 -1+2 v 


02 = 


(1 + 2v)B± |F| 1+2,/ 


+ • • • + (|y|"- v + + 


2(3 - 2 1/) F 2 

+£ 2 /i* +1/ B| : 52±|F|^ + ‘ / + «y-'*i* a± |F|H' + • • • , 
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^eV~0 - tV" ^ 




2i/(3 - 2*/) F 2 


4fc 


ys 




4i/ 


( 1 - 2 */) 


(1 — 2i/) 

|y|* -1 ' + 


a-" 


(3.15a, 6) 


frfo _ 2 -l+2i/ 

dr ~ V 


(i + 2tQ a|B+[ 2 iri 1+2,/ 


dr F 2 


+ 


dr 


= eV 1+2l/ .0 - e 2 /! 2 " Jl a|5±|2 


2i/ dr y 2 


+c 2 «// _4+21 ' 


q |y|l+2*'l 

|(1 - 4^ 2 )(1 - 2^(77 - l)[^ ± pL^__ 


+ 


(3.16a, 6) 


, 3 -|+3t/(l + 2i/)(l — 2i/) 

* = - £ '‘ 1 2(3 + 2v) 2 |K| 


Bi 


f — 3i/ 


+ 


<f>3 = t 3 fi~^ +3l/ .0 + • • • , 

(3.17a, b) 

Note that above, for the sake of brevity, we have only included the more important and 
illustrative terms with regard to our aim of deriving amplitude-equations for B. It should be 
noted that the above asymptotes, in general, only contain the leading term (for Y — ► oo) at each 
order. Further, it should also be noted that these (incomplete) asymptotes may, depending on 
the relative sizes of e, k and fi, be disordered as written; it is sensible to postpone a discussion 
concerning the relative sizes of these quantities until §5, when the sizes of the competing nonlin- 
earities are known. Here d±,c± and q± are constants, as introduced by Churilov & Shukhman 
(1988); the first two are in fact expansions, each term of an order yet to be determined. The 
constants <7 2 ± and /i 2 ± are defined in Appendix A of this paper. 

It is worthwhile to consider these asymptotes a little further. We see that the leading terms 
in the ^-expansions are all zero; in fact, such behavior is in full agreement with the theorem 
of Miles (1961). In the next section, we shall see that this fact results in the ‘biggest’ cubic- 
nonlinearity not contributing to the amplitude equation; the same result is also found when J = 
1/4 (y = 0) and is responsible for the quintic nonlinearities occurring in the evolution equations 
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derived by Brown & Stewartson (1978), Churilov & Shukhman (1988) and in this paper. As 
pointed out by these previous authors, for a ‘contributing’ jump, it is in general necessary for the 
(^-asymptotes to be non-zero. They note that this condition is satisfied if the complementary- 
function terms in the <f> 2 asymptote are non-zero (thus resulting in a ‘contributing’ quintic 
nonlinearity). Churilov & Shukhman (1988) also note that if viscosity effects are included in 
in the original formulation of the problem, then the above condition is satisfied as long as the 
Prandtl number rj ^ 1 (thus resulting in a ‘contributing’ cubic nonlinearity). In addition to 
these two so-called symmetry breaking mechanisms, it is clearly obvious from the third term 
of (3.14b) that we have an additional mechanism for non-zero u, namely non-zero J\. This 
novel mechanism is made possible by the two Frobenius roots, 1/2 ± \v\, not differing by an 
integer. As far as the authors are aware, such a nonlinear, non-equilibrium critical-layer has 
not been the subject of any previous studies. This ‘new’ mechanism is essentially responsible 
for the ‘non-viscous’ cubic jump term in our later amplitude equation when v is negative (later 
referred to as the ‘ Ji-cubic’ term). 

The fourth symmetry-breaking mechanism possible here corresponds to that cleverly identi- 
fied by Churilov & Shukhman (1987) in the unpublished Appendix B of their paper. Ironically, 
there is no mention of it in the sequel paper (Churilov Shukhman, 1988); however, it is of 
lesser importance for disturbances of marginally unstable flows (see later for an extended dis- 
cussion on where each of the four possible nonlinearities in the evolution equation for B are the 
most significant). The important thing to note at this stage is that the term involving h 2 ± in 
the <f > 2 asymptote is not at one of the orders to be considered, in §4, for the preceding term in- 
volving c 2 ±. In the next section, we shall see that this symmetry-breaking mechanism provides 
the ‘non-viscous’ cubic jump when v is positive (later referred to as the outer-complementary- 
function [OCF] cubic term). 


4 The Critical Layer 


As usual in such nonlinear studies, the main pu rpose of this section is to calculate the 
second relations between &[+ and b^} (the first being given by the solvability condition, (3.9b)) 
and thereby obtain the desired nonlinear evolution equation(s) for the disturbance amplitude 
B(t). Following Brown & Stewartson (1978) and Churilov & Shukhman (1988), we define new 
functions 'P, P and $ where 


^=L 2 y 2 + iIr, p = p 0 -pY -2/x" 1 iV, $ = tf-(l+2i/)P. (4.1a - c) 

Inserting these into the governing equations (2.1), (2.2) gives 

= 5(1 - 2i')$r - H~ 2 {Vy, yy - ^ - *,) + r,Kp- 3 ^YYY, (4.2 a) 

%-.„)* = -/P 2 {*k, *}* - - *«) + W~*$YYY + (17 - l)K/i- 3 *yyy, (4.26) 
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where {a, b}* = a x by — ay 6*, the operator 



and it is assumed that /c/z -3 < 1 i.e. viscous effects are not large enough to affect the operator 
N x at leading order. 

In the rest of this section we proceed to solve equations (4.2) for the relevant lower order 
terms of the lower harmonics. In fact, we need to consider the cases v > 0 and v < 0 separately; 
this is because the largest ‘non-viscous-in-origin’ cubic term in the corresponding evolution 
equations is different in these cases. When v < 0, this term is proportional to J\\ but when 
v > 0, this term is related to the complementary function of the second harmonic of the inviscid, 
outer problem. Note that when v > 0, the second and third terms of the V>i asymptote (3.14a) 
art ordered and thus the cubic nonlinearity proportional to J i would be insignificant compared 
to that proportional to Note that the other two significant nonlinearities (a cubic, due to 
viscosity, and a quintic, due to complementary function terms of the critical layer solution for 
the second harmonic) are possible in either case. 


(i) The case v < 0 

The solution is again constructed in the form of a Fourier series in x; we expand the 
fundamental, zeroth, second and third harmonics, respectively, as follows 

1 = e / x2+‘''f ( 1 1) + • • • + e 3 /z-3 +3,/ $i 2) + • ■ • + e^ + ‘ / ’P ( 1 3a) + • • • + + • • • 

+c 3^-|+3^( 4) + . . . + + • • ■ + e 3 K/z-^ +3,, tf S 6) + • ■ • + eV^ 5 ^! 71 + • • • , 

+ . . . + eV 2 ^S 2) + ■ • • + c 2 K/z- 4+2 ^i 3) + • ■ • + + ■ • • , 

+ • • ■ + ’ • + e 2 K f i- 4+2l '¥^ ) + ■■■ e 4 fi- 4+4t/ ¥j ) + • ■ • , 


(4.4a — d ) 

and similarly for the $,’s. These expansions are not necessarily completely ordered (depending 
on the sizes of e, k and /z); moreover, only those terms crucial to deriving the evolution equations 
have been included. The scalings follow directly from the outer asymptotes (3.14-17) and/or 
by considering the process of harmonic generation. 
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4.1 0 (e/z 2 +,/ ) of the fundamental 


Upon writing 


Nl x = 


' d_ 
<dr 


+ iklY 


) 


d 

dY 


iklx , 


we find that at this order equations (4.2a,b) give 

JV., ,}+„>* i" = 5(1 - 2 »)M?\ JV-Ml--)*” = 0 ■ 

A solution of these equations which matches to the outer solution (3.5) is 


y[ l) = W(t, Y) = 


(1 + 2v)T(\ + v)u_ 3l/ 

47rA;2 + ' / 

*i 0 = 0, 


L 


t-%- v B{T - t)e~ ikYt dt, 


(4.5) 


(4.6a, 6) 


(4.7a) 

(4.76) 


where the contour C is shown in Figure 4.1a. The function W(t,Y ) has a single asymptotic 
representation in the lower half-plane (— ir < argU < 0) 


W{t,Y) = B{t)Y^ u + 0{Y~^ v \ as |y|->oo. 


(4.8) 


Later we derive evolution equations for the amplitude B{t ) but for the moment we can regard 
it as an arbitrary function that satisfies the requirement B{t) — ► 0 as r — ♦ — 00 . Matching with 
the ‘outer’ asymptote (3.14a) yields 

B _ = r 1 ~ 2l 'B + ( B + = B); ( 4 . 9 ) 

this result has already been referred to and used in the previous sections. 


4.2 0(e 2 ^ _1+2 ") of the zeroth harmonic 


At this order equations (4.2a,b) yield 

= ik(WW - WW')\ = 0, (4.10a, b) 

or or 

where the overbar again denotes the complex conjugate and prime denotes differentiation with 
respect to Y. A solution of the above which matches to the outer solutions is 

*?’ = *S ,) = 0. (4.11a, b) 


12 



4.3 0(e 2 /i -1+2l/ ) of the second harmonic 


At this order equations (4.2a,b) yield 

^ ( M®a 1)=ifc (l”2^)®f ) + * fc ( , ^-^ a ) and N Hk-^ ] = °- ( 4 * 12a > 6 ) 

We note that, as the right-hand sides of (4.12a,b) do not involve W, ^ ^ an( ^ ^2 ^ have 
unique asymptotic representations as \Y\ -* 00 (in the lower half plane of complex F) of the 
form 


vpW „ ± + ( m^F* - " + n i 2 ) Yi +,/ ) + • • • , $ 2 1} 


4 vm^ T ,i_ 


(1 — 2u) 


Y*- y + 


(4.13a, 6) 


3 ( 1 ^ 

Matching with the asymptotes (3.15a,b) requires that c± = e 2 /x~ 2 +3l, ci±, so that ci + = m 2 , 
Cl _ = n 2,/ c 1+ q + = n 2 1} and pFcx-q- = t'- 1 “ 2, 'n£ 1) . These lead to the relation 

C\-(q- - i~ 4l/ q+) = 0. However, as the adopted model flow is symmetric, q + = q- and so, for 
v / 0, ci_ = 0 = m 2 ^ = Thus 


$i 1} = 0 


(4.14a) 


and the solution to (4.12a) is 


= 


i 1 ~ 6t 'k 1 ~ 2u (l + 2i/) 2 r 2 (| + t/) 


32tt 2 


f dti f dt 2 f dtB(T—t—ti)x 
Jc Jc Jo 


5 ( r _t_t 2 )(t 1 t 2 )-|-*'(i 1 -t 2 ) 2 (t 1 +t 2 )^(2t+t 1 +< 2 )-t- ,/ e-’ w ( 2<+<1+<s ). (4.146) 


4.4 0(e z fi~^ +3l/ ) of the fundamental 

From (4.2a,b), 

JV, (1+H *W = ifc-(l - + 2 |I‘V - *?> V - #?>>), 

* 2 ^ 

N ui _^? = 0. (4.15o, ft) 

Due to the earlier $ terms being zero, the right-hand side of (4.15b) is also zero. Thus 
e^V, where is a constant and V satisfies N- l ^i_ l/ ^V = 0 with the asymptotic representation 
V ~ Y\~ v as |F| -v 00 . Any non-zero choice of e^ 2) in (4.15a) results in a Fa - " term in the 
asymptotic representation of ^ and this contradicts the theorem of Miles (1961), from which 
we know that the expansion of tf> a as y -4 0 cannot involve a term proportional to y%~ v . Thus 
it follows that $ 1 ^ = 0. 
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We also note that the right-hand side of equation (4.15a) contains W and so \&i 2 ^ does 
not have a unique asymptotic representation as |T| — >■ oo. However, because the asymptotic 
representation of ' does not contain a term proportional to Y J - ", it is impossible to obtain 
non-zero at this order. Thus the strongest nonlinearity does not contribute to the evolution 
equation. 

The solution for \ ’ is needed for later calculations regarding the quintic nonlinearity; the 
solution of (4.15a) is 




2 -|-7«/£§-3«/(i _|_ 2i/) 2 T 3 (| + v) 
32tt 3 


f°°dt [ dh ( dt 2 f dts 
Jo Jc Jc Jc 


B(t - t - ti)B(r - t - t 2 )B(r - t - t 3 )(M2*3) * u {t 3 -t\- -t-ti- t 2 )-%~ v x 

^3^l(h — ^3) 2 + ^2(^2 — ^3) 2 + ^ 2 ( 2^3 ~ tl ~ ^2)| _ t\^ 2v t 2 {ti + t 2 )~*~ V X 

(h-t 2 )i~ , '(t l +t 2 + t3)(2t3-ti -t 2 )F^---v,~- 1 -tl) 

(4.16) 

Here H{x) is the so-called Heaviside function: 



1 : x > 0 
0 : x < 0, 


and F(a, b ; c; z) is the hypergeometric function (see Erdelyi, 1953; Abramowitz & Stegun, 1964). 


4.5 0(e/J+") of the fundamental 


A major difference of our problem for general 1 /, compared to the u = 0 case considered by 
Churilov & Shukhman (1988), is that there is no linear contribution to the evolution equation 
from the critical layer. This can be seen from the analysis of this and the next subsection. 

At order t^ v terms proportional to J x first enter; equations (4.2a,b) yield 


^,a + ,,*i 3a) = 5(1 - 2 K)i»r - 


(3a) 


2ikJ\ 


(1 + 2u) 


W, 


> 0 “) - 2iUi 


which have solutions 


N (TT^y 


*f ) — ^HMn0.y). * ( i a) = -^T^)- 


(4.17a, b) 


(4.18a, b) 


Note that for definiteness, we choose to include the complementary function alongside the 
nonlinear term which is going to be balanced with at this order. 
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4.6 0 (e/j 2 _ ") of the fundamental 

Here we find that (4.2a, b) yield 

= °- 

These have solution 

= *f 6) = 0, (4.20a, 6) 

where once again we shall include the associated complementary function alongside the non- 
linear term which is going to be balanced with at this order. The solutions (4.20) imply that 
the whole of the linear contribution to the evolution equation comes from outside of the crit- 
ical layer; note that the logarithm occurs in (4.18a) and not in (4.20a) i.e. any jump induced 
by the presence of the logarithm does not occur at the desired order, 0 (e/i* - "), necessary to 
affect the relationship between b[ 2 } and 6$. The authors have verified that this is in fact the 
case by computing 4 near —neutral’ linear growth rates from a numerical solution of the Taylor 
Goldstein equation (2.2) and comparing them with those predicted analytically by equating the 
right-hand side of expression (3.9b) to zero. 


4.7 O (e 2 /i 2l 0 of the zeroth harmonic 

At this order we only need to calculate 3>o^i (4.2b) gives 


with solution 


— = , tUl AWW - WW')\ 

dr 0 K1+ 2^) 

(4.21) 

,(2) 2Jl IW'I 2 — Jl 

0 ” 1/(1 + 2uf 1 1 ~ 1/(1 +2^) 

(4.22) 


4.8 O (e 2 fi 2u ) of the second harmonic 


Equations (4.2a,b) yield 

N 2>c _ +l/) ¥ 2) = t*(l - 2i/)'*? ) + ik{W"V S 3a) - 2 w'¥? aY + w^S 3a) ") - n 4 +^ 2 1} , 


(1 + 2u) 


"mm*?* = -W^r w " - ^ - <rf (4 - 23 “' 6) 

For further calculations only $ 2 2 ^ is needed; using the identity AT 2> (i +1/ ) = N 2 ^i_^ — 4 ikv and 
comparing the right-hand sides of (4.12a) and (4.23b) it is easy to deduce that 


$(*> = h — ^ . 

2 i/(l + 2v) 2 


(4.24) 
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4.9 0(e 3 /z - § +3t/ ) of the fundamental 


At this order the cubic-jump proportional to Ji emerges. The governing equations yield in 
concise form 


N, 


<■ yTf( 4 ) — In 2ikJi y(2) p(4) 

l,(i+ 0 W l - M)** 9 ! Mj. 9 .A* 1 +Kl ’ 


2 V ‘ (1 + 2u) 

2ikJi 


Hui-.,*? = + ^ 4) ’ ( 4 - 25 ‘- ») 

where, in particular, 

4 4) = + 2?1 8b) "®? ) - *<? a) W + *?>>- ^2 2) "^ - 2$5 2) 'TF') 




(-2W*P + - W9& 1)U + ^2 Y 'w + * { } Y W ) 


(4.26) 


“ i/(l + 2 v) 

To derive an evolution equation the asymptotic representation of VP^ is needed. This is 

$< 4) = c£ ) Y> +u + D < £ ) Y±- v + 0{Y~* +U ) as Y — ► ±00 , (4.27) 

where, in particular, 


D ? - D - = -sfe L dY r dt - *• y)i 


- ikYt 


(4.28) 


Note that the term proportional to in (4.25b) only contributes to the jump — Cl 
After substituting for and some calculation, we find that 

r>(4) r>( 4 ) * 1_4l/ A; 3-4l/ J \ ( 1 + 2i/)T 3 ( | + i/) t r t . 

D *- D -= 32*W(1 - V) J c dh J c dh J c dt 3 (tM 3 ) > x 

^ 5 (r + <1 — tz)B(r -f- <2 ~ tz)B{j + t\ + <2 — 2 tz)t 2 tz{tz ~ ^2)^3 4 - t\ ~ ^2)* 

(^3 ~ t\ ~ t 2 ) 2 v H{tz — t\ — t 2 ) 

-f-- 1 f dtB(r + ti + t — ^3 )B(t + t 2 4 t — tz)B(r + t\ + ^2 T 2 t — 2^3) x 

4 Jo 


.(4) 


(ti - f 2 ) 2 (*i + h)^ v (2t + <! + h)-*~ u {2h -2t-ti- < 2 )x 

> 

(2t + ^ 1+^2 + ta)(t 3 — 2 t — ti — t 2 ) 2 v H(tz — 2t — ti — t 2 ) 


(4.29) 


Matching the ‘inner’ asymptote (4.27) with the ‘outer’ asymptote (3.14a) leads to the rela- 
tions 

z4 4) = &$£+, £>L 4) = z 1 - 2 ^! 2 .^-, 

which combine to give 

£>(4) _ £,(4) _ B ( b (2) _ *-4^(2)^ 5 (4.30) 

where the subscript J x denotes the cubic part of the total jump which is proportional to Jj. 
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4.10 0(e 3 /x"* +3 *') of the third harmonic 

At this order the governing equations yield 

jV 3 ,(i + „)*i 1) = ^(1 + 2 1 /)$| ,) + ik (W'w - 3*?*' w' + 2¥< I W). 

)*?’ = 0- (4.31a, 6) 

The right-hand sides do not contain W V and so have a unique asymptotic representation. Match- 
ing with the solution outside the critical layer gives $ 3 ^ = 0, whilst solving (4.30a) gives 


$ 


(i) 


tf~ lll 'A:f~ 3 ‘ / (l + 2i/) 2 T 2 (| + v) f 




dt dt a dt 2 dt 3 B(r -t-tjx 


B{T-t-t 2 )B{T-t-t 3 )t\ +v t;* v t^ , '{t 1 + t 2 )-2-‘ / {h-t 2 )2 u x 
(h +t 2 + f 3 )* + "(3i + h+t 2 + i 8 )“ l "*'(h + U- * 3 )(h + h ~ 2< 3) x 

f(-\ - -'-i - f ’ ! - b *ly‘ tY{M+M,M) H(h - k). ( 4 . 31 c) 


4.11 O (eV -4+4l/ ) of the zeroth harmonic 


At this order it is only necessary to determine Equation (4.2b) gives 

JUW = 0, (4.32) 

which has the solution = f(Y). However, choosing a non-zero f(Y ) is merely equivalent 
to taking a different density profile for the unperturbed flow and thus, without any loss of 
generality, we choose to take f(Y) = 0 and = 0 . 


4.12 O (e 4 /LT 4+4l/ ) of the second harmonic 

At this order the equations (4.2a, b) can be written in the form 

A 2i (£ +t/ )'P 2^ — *&(1 — 2 i/)$ 2 7 ^ = ikRff = ik(Rn + R02 + #31)1 

N 2 ,(i- v) *P = 0, (4.33a, 6 ) 
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where 


R n = (W'l'f 0 " - 2 W'qW + W "<p[ 2) ), Rm = 2(V ( J ) "¥ 2 ) - >') 

and R 31 = (3TF"^S 1} - 2lF' r - W^")- (4.34a - c) 

Following Churilov & Shukhman (1988), is written as the sum ’J’P = 'I'jtv + ® 2 l °f 
a particular solution (-^.(i+iO^/v = ikRjf) and a general solution, \& 2 L? °f the associated 
homogeneous equation. Defining 

we find that 

$< 7 > = /i - 2 l/ co 0 (l /2 - t/, 2 k) [ rl +v C(r - t)e~ 2iktY dt , 

= coo(l/2 + i/,2Jfc) / r%- v D{T - t)e~ 2ikiY dt 

Jc 

—fi~ 2v — Coo(l/2 - t/, 2Jfe) J t~*+ v C(r - t)e~ 2ikiY dt, (4.35a, 6) 

where C(r) and D(r) are arbitrary functions which tend to zero as r — *• — oo, and the fi~ 2l/ 
factors have been added for later convenience. As Y — ► ±oo, 

= M±Y* +V + 0(F-2+ 1 '), 

where 

= oT^T7) L iY S~ iiR ” {T *°- (436) 

Thus the asymptotic representations of * and $2 ; have the forms 

^ 7) = [D(t) + M±(r)]K^ _ ^1c{t)Y*- v + • • ■ , 

$5 7) = pT 2v C{t)Y %- v + • • '• . 

(4.37a, b ) 

To match with the solution outside the critical layer, in (3.15) it is necessary to set c± = 
e 4 ^ - l +3l, c 2 ±, leading to the relations 

c 2+ = i 1 ~ 2l, c 2 - — C(r), c 2 +q+ = D(t) + M+(t ) and c 2 _g_ = i 1+2l/ [D(T ) + AT_(t)]. 

Upon setting q+ = q = q-, we find tha,t 

C'(t) = ~ l ~^~ csc (■ y)( M + “ M -)> ( 4 -3 8 ) 

and, as M + — M_ ^ 0 , this fully determines 7 ^ 0 . Thus, the symmetry $ = 0 has been 
broken and a ‘jump’ will occur at the next order in which $ 2 ^ appears in the right-hand side 
of the equations. The explicit form of M+ — M_ is given in Appendix B of this paper. 
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4.13 0(e 5 /i _ ^ +5,/ ) of the fundamental 


At this order, the governing equations are 


Wum** 7 ’ = ?*< + 

JV 1 , (H) *< 7) = -iV (*P" *?> + 2*‘ T) '*i ,) ') = 

The asymptotic representation of has the form 

$! 7) = C { l ) Y^ v + D ( ± ) Y>- 1/ + 0(r-* + ‘'), 

where, in particular, 


D i 7) - Z><. 7) = 


i%- v k*- v f°° IT , /°° 


1°° dY rd„i-^\r-t,Y)e-‘ 

k — V) J- cc */0 


-»*yt 


2^(1 - n) 


From (4.7a) and (4.35a), 


(4.39a, b) 


(4.40) 


= J-^(l-4^)cos(^) / , _ (j)e _ il y (a ,- b)c , (T _ tl)B(r _ y, 

2l * Jc JC (4.41) 

and simple manipulations give 

^+ 7) - = C *J 0 t '~ vc ( T - *) 5 ( r ~ 2i ) dt ' c * = „)r(i 4- |)rd - 1) 


2ii -3, 'fct- , '(l - 2 i/)ttI „ n 

. . r'r/l _ t 4 ' 420 ' 6 ) 


The quintic contribution to bQ is found by matching the above asymptote with the outside 
expansion of the fundamental near the critical layer. This process yields 

- B(6' ! + > - rn<?)„ (4.43) 

where the subscript 5 denotes the ‘quintic’ part of the total jump. From (4.38), (4.42) and (4.43) 
we obtain 

B( 6*2 - = -&scff) r t [M + ( r - 1 ) - M.(r - ()] B(r - 2t)dt. (4.44) 

/(7 Z */0 
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4.14 The cubic nonlinearity due to viscosity 


The remaining terms in the expansions (4.4), proportional to k, are driven by dissipative 
(viscous) effects. They are not of major concern in this paper but are included for completeness; 
the equations which they satisfy, their solutions and further analysis can be found in Dando 
(1993) (see also Churilov & Shukhman, 1988). Providing the Prandtl number is not unity, a 
non-zero jump, — D ^ — Z?L 6 ^]o + [D+^ — dL 6 ^] 2 say, arises at 0(t 3 Kfi~ 1 £ +3i/ ) where 


[ z 4 6) -£- 6) ] 0 = 


■1—4^5— 4^ _ 1)T3(1 + „)(! + 

32i/7t 2 (3 — 2v)T{\ — v) 


f dt\ f dt 2 j dt 3 

J C J C J c 


B(t — t\ + ^2) (^1 ^2^3) 2 ^(ti — t 2 — £3) 2 ^1 "■ ^2 ) 2 1^2 j dt^B{r — 2t\ + <2 + ^3 ^4) x 

B(t - ti + <3 - t 4 ) ^(1 + 2v){t \ + <2) - 2(3 - 2v)tit 2 {t 2 - h) 2 ^ 

—B(t — 2t\ + t 2 + t 3 )B(r — t\ + ^3)^1 t 2 ^(ti — t 2 ) 3 ~ 2< 3 (<i — t 2 ) 2 + H{i\ — t 2 — t 3 ), 
and 


[D? 


( 6 ) 


■D% = 


•1-4, A .5-4,( t? _ i) r 3(I + „)(! + 2 V f 

256^(3 - 2v)Y{\ - u) 


f dt\ f dt2 f dt% f dt 
Jc Jc Jc Jo 


B(t — 2ti + tz + ^3 + 2 t)B{r — t\ 4- £3 + i)B{r — t\ + ^2 -f 0(^i^3) 2 v x 
(^1 — ^2 “ ^3 “ 2^)2 ^(^3 — ^2) 2 (^2 + ^3) 2+i/ (2f + ^2 “H ^3) 2 ^(2^1 — i>2 ^3 “ 2t)x 

^2tf + (2f + 1 2 + t3) 1+2l/ (t 2 — t 3 ) 2 (t 2 + tz) 1 2v + (2t + t 2 + t3) 4 ^/f(ti — t 2 — t 3 — 2 1). (4.45a, 6) 

Matching with the outer inviscid solution requires that 

Df - Z)l 6) = B{b?l - i~ 4l/ b ( Y-)v, (4.46) 


where the subscript v denotes the cubic part of the total jump which is due to dissipative effects. 

It is worthwhile to consider the effects of viscosity a little further at this point. At the start 
of this section, it was pointed out that we assume « <C p 3 , so that the effects of viscosity do 
not enter the crucial critical-layer operator N x at leading order. Instead, the effects of viscosity 
occur as inhomogeneities at certain lower orders of the hierarchy of critical-layer equations. 
This approach is entirely rational as long as «//i 3 <C 1. However, to derive amplitude equations 
valid over a larger range of k values, one must follow the approach, introduced by Haberman 
(1972), of introducing a new parameter, A h say, where 

A# = /c//i 3 , (4.47) 
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is taken to be order one during the critical-layer analysis. The operator N x would need modi- 
fying: 

N * — * N x - (4.48) 

and hence all the critical-layer solutions would be modified. We have chosen not to adopt this 
approach here as it would complicate the presentation of our theory and analysis; whilst it is 
not really necessary so long as A// is reasonably small (see Goldstein & Leib, 1989). 


(ii) The case v > 0 

Clearly the viscous-cubic and the quintic nonlinearities are still possible when v is positive. 
Moreover, the analysis and the evolution equations are exactly the same as for the negative v 
case. However, when v is positive e//5 -v > e^ JrV and the ordering of terms in expansion (4.4a) 
strictly needs to be changed so that the term is lower order than the term , It is now 
no longer rational to balance the cubic term formed from the a term, with the ; term. 
Thus, the J \ -cubic nonlinearity considered for negative v is no longer a possibility. 

There are two other possibilities that we can consider for positive v. Firstly there is the 
cubic formed by the * term, balancing with the term. However, this among other 
things, necessarily leads to the fully nonlinear critical layer problem. Alternately, there is 
another cubic nonlinearity ( referred to here as the outer-complementary-function [OCF] cubic 
term) similar to that considered by Churilov & Shukhman (1987) (for unity Prandtl number) 
in the unpublished Appendix B of their paper. This cubic nonlinearity was not considered for 
negative v as the Ji-cubic nonlinearity is always larger. It arises from considering the part of 
the complementary function term in the asymptote of the second harmonic outside the critical 
layer (equation (3.15a)) to be at an order fixed by the outside, inviscid problem; rather than 
just considering it at orders fixed by the process of harmonic generation inside the critical layer. 

It is necessary to consider two additional terms in the (slightly re-ordered) expansions 
(4.4), namely terms e 2 $ 2 ^) i n th e expansions of the second harmonic and terms 

e 3 /x -1 (^[ 8 \ $ 1 ^) in the expansion of the fundamental. 


4.15 0 (e 2 /i 2 _l/ ) of the second harmonic 

At this order equations (4.2a,b) yield 

N u i +1/) ^ 8) = (1 - 2 v)ik$ { ? ] + ik(WW" - W' 2 ), and W 2t( i_,, ) *? ) = 0. (4.49a - b ) 

The solution for i® written 

= coo(l/2 - 1 /, 2k) [ r f +,/ n 2 (r - t)e~ 2iktY dt, (4.50) 

J c 
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where the function n 2 (r) is determined below. For large \Y\, 

$< 8) ~ n 2 {T)Y*- v , ~ m 2 {r)Y^ v - - ~^ -n 2 (r)F?-", (4.51) 

where m 2 (T) is a suitable arbitrary function. Matching with the expansions outside of the 
critical layer yields 

< 72 -=*' 1+2 ^ 2 +, h 2 .=i 1+6, 'h 2+ and B 2 h 2 + = - 

From the relations (4.52a-c) and (A.5d), it follows that 


, . 1 


n2 ^ T ' (1 - 2v)a a sin(7 tv) 

and thus as is non-zero, it will lead to a ‘contributing’ nonlinear jump. 


— 2v) 

4j/ ” 2 ‘ 

(4.52a - c) 

-S 2 (r), 

(4.53) 


4.16 0(e 3 /z _1 ) of the fundamental 


The governing equations at this order are 

iV li( i +I/) tf ! 8) = ^(1 - 2v)ik$ ( ? ) + ik{2W'v? - W - WV ( 2 8) "), 

N, = -ik(W^ Y ' + 2W 1 $f Y ) = - ikR f\ 

The asymptotic representation of as Y — ► ±oo has the form 

vlrj 8) = C^Y**" + D < £ ) Y5-‘ / + 0{Y~^ V ), 


(4.54a, 6) 


where, in particular, D+ 8) - D_ 8) is given by the right-hand side of equation (4.40) but with 
replacing R 2 \ Simple manipulation gives 


= _2t*"ik(l - iv^costjv) t dh [ dh g {T _ (])nj(r _ h)e -^-u)Y t -h- t 
2 7T Jc Jc 


t 2 
C 2 


(4.55) 


and we find that 

1 ,( 8 ) n<*> 64»5~ H 'A:!~‘'( 1 - 2i/)r(^ + i-)irl /”* ti-'BjT - 2t)ni(r - t)dt (4.56) 

D * ~ - " »r(} - |)r(J + 1) h 

Matching with the outer, inviscid solution requires that 

Df - D ?> = B(6«,5 - (4-57) 

where the subscript oc/ denotes the cubic part of the total jump which is due to the outer- 
complementary-function terms. 
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5 The Evolution Equations 


The four possible forms of the evolution equation for B{t ) simply correspond to matching 
the solvability condition (3.9b) with each of the four jump expressions (4.29), (4.43), (4.45) and 
(4.56) in turn; for the sake of succinctness, we do not explicitly quote these evolution equations. 
Instead, it is more beneficial to the discussion here to form one ‘composite’ evolution equation 


2il\ 


(i+i 4 n 


dB Jihk 

«7 + ir l “ WB 


= - c lvB(b { il-i- 4,/ b { ?]) 


= -2 - i~ 4v b^l)j i H{—v) + eV l+ "(&i+ “ i~ 4u b ( yl)oc}H{v) 

+^ f r 7 * i -K(b^ - m<l>), + £ v-w(t«5 - mil 1 ),}, (5.i) 


where the explicit forms of the nonlinear quantities appearing on the right-hand side of this 
equation have been derived, And are quoted, in §4; again H denotes the Heaviside function. As 
the composite equation (5.1) contains the sum of the four (individual) nonlinearities, it has the 
advantage of being valid for all values of k, e and /i that lead to an unsteady (weakly nonlinear 
and weakly viscous) critical-layer. 

In this section we shall discuss the parameter ranges of validity of each of the four possible 
forms of the individual ‘base’ evolution equations for B(t), before presenting some numerical 
results for the ‘ J x ’-nonlinearity case. This section is concluded with a discussion of the expected 
solution properties of the other (base) evolution equations. Further conclusions are drawn in 
the next section. 


5.1 Parameter ranges of validity 


The range of validity/application of each of the possible base evolution equations is sum- 
marised in Figures 5.1a,b, for v < 0 and v > 0 respectively; the governing balances determining 
the dashed curves separating the three parts of region III follow immediately by comparing the 
sizes of the four terms on the right-hand side of equation (5.1). As an example, we note that 
(when v < 0) the nonlinearity proportional to J\ is equal in size or larger than the quintic 
nonlinearity if eV 3+4 ‘' > eV _7+6 " i- e - e ^ ^ 2_ " • We see that in each case, depending on the 
relative sizes of /c,c and //, three of the four base evolution equations are applicable. Also plot- 
ted are the expected evolutionary paths of the disturbance for each case; this shall be discussed 
in more detail later but at present it suffices to note that, based on our assumptions concerning 
viscous-spreading effects resulting in an unstable linear disturbance mode approaching a later 
neutral state, initially our disturbance will lie in the bottom right-hand corner of region Illb 
in each case. This partly justifies our relegation of viscosity to lower order effects, as well as 
indicating that the base evolution equations with the J x -cubic and the OCF-cubic nonlineari- 
ties deserve the first attention. In the next subsection we present numerical calculations for the 
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former case ( 1 / < 0); the evolution equation with OCF-cubic nonlinearity {v > 0) is discussed 
in the final subsection. 


5.2 Numerical results for the Ji-cubic 


Recall that if, and only if, v < 0 then a cubic nonlinearity proportional to J\ is possi- 
ble. In fact it is the largest nonlinearity in the (composite) evolution equation when p k* 
and e <C p 2 ~ v \ thus, as mentioned in the previous section, it corresponds to the first nonlin- 
ear evolutionary stage for an unstable linear disturbance whose growth rate is diminished by 
viscous— spreading effects. In this subsection we present some numerical results for this case. 

To ease numerical calculations the jump expression (4.29) is transformed into kernel form, 


cW-Di 4 ' = ±M [“iss 1 - 4 - ['dv<,'- 2 *G(c)B(T- s )B(T-a S )B(T-(U-<7)s) 

x 4 (j — v) Jo Jo 


- B(b$ - i-n 


(5.2a) 


where 


G{a) = 


r 2 (| - v)0 + <r)“* 


— i— 1/ 


(i + 2,)r ( i - 2.) { 2(1 + ^ - 2v)F 4 (I - y ’ 1 + ^ I + ^ 1 ~ 21,5 (TTg)) 

-(1 + <J + <t 2 )(3 - 2u)F\ \ + v \ + ^ 2 - 2v 'y *’> (TTv)) 

cr(l + <r)(5 — 2i/)(3 — 2v) (1 3, 1, 0 0 <7 \\ 

+ ~ fco F '(-2- l/ '2 +I/ '2 +V '*- 2 v ^0T7))] 

+(1 + <r)(l — (t) _2i/ x 

r \ *-§-"(1 - OWl - at )~%+ Zv , / 1 1 v 3 v 2 ,\' 

l dt- J- 4 (2-cr+<?t-2cr t 2 ) . 

Jo (1 -f- crt) 2 "*’ l/ (l — <j 2 i)*+ v V 2 4 2 4 2 / 

(5.26) 

Here F\(a , 6, c, d\ x, y) is the hypergeometric function of two variables (see Erdelyi, 1953; Churilov 
k Shukhman, 1988). The evolution equation for this case can then be written in the form 

= llB + 72 r ds s 2 - 41 ' Pda a l - 2v d{<T)B{T - s)B(t - (ts)B(t - (1 + a)s), (5.3) 

OT JO Jo 


where 


Jik^i t\1 + 2u\ a M , 0 ,.^ 

7i = 72 = 77-77175 / \ t and G(a) = -sgn(l + 2t/)G. 


2I\ 


8t/r 4 (| - u) cos(7ri/)/x 


(5.4a — c) 
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Let us consider the range -2 < v < 0; in Figures 5.2a, b -71M and -/faM are plotted 
versus v for this range. These figures illustrate that both 71 and 72 are positive in this range 
(recall that, since we are considering unstable modes, J\ < 0). The lower value v — —2 chosen 
here for this numerical investigation has no special significance, it was principally chosen so 
that we calculate numerical results for a complete period of the trigonometric functions present 
in the analytical expressions. Note though that, for their study of the related neutral inviscid 
Gortler modes in a three-dimensional boundary layer, Blackaby Sz Choudhari (1993) find that 
this range of ^-values is appropriate. Moreover, it may also be argued that for increasingly 
negative values of v (corresponding to the Richardson number becoming increasingly negative) 
the now destabilising effect of stratification means that the growth rates of the more unstable 
linear inviscid disturbances are just too large for viscous-spreading effects to damp significantly , 
thus rendering the weakly— nonlinear theory inapplicable. 

It is convenient for numerical calculations to introduce a so-called ‘logarithmic time (see 


Churilov h Shukhman, 

1988; Shukhman, 1991) 




(5.5a) 

having set 

f 

Q? 

II 

k 

(5.56) 


where the constant B 0 is chosen such that 6(r) 1 as r -► -00. This is done in order to 

reduce the number of parameters in the evolution equation. Equation (5.2) now has the form 


— = [' daK(a) I" dxx 2 ~ 4v e~ x b(Te~ x ^ 1+ ^)b(Te~ ax ^ 1+ ^)b(Te~ x ), (5.6 a) 

dT Jo Jo 

with initial condition 

b(T = 0 ) = 1 ; ( 5 - 66 ) 

where the kernel 12(/ 

KM = (1 + ^ gW- (5 ' 6c) 

In Figure 5.3 K(a) = K(ar\v) is plotted, versus cr, for a few representative v values. It is 
interesting to note that K(a) is always negative (0 < cr < 1); this appears to be the case for all 
—2 < v < 0; cf. the viscous-jump kernel for the v = 0 ( J = 1/4) case considered by Churilov & 
Shukhman (1988) which changes sign. Thus it is possible to deduce all the qualitative results 
of the solution properties of the evolution equation (5.6) from results plotted in this figure (see 
the discussion in §5.2). However, we still chose to solve the integro-differential equation (5.6) 
numerically for completeness, as well as to obtain actual quantitative results. 

In Figure 5.4, we present the results of a numerical solution of the evolution equation (5.6) 
for two representative v values. The results show that 6(T) oscillates with a fast rising amplitude 
as T increases. The period of these oscillations is seen to depend strongly on u: for v values 
corresponding to the smaller x 2 ~ 4l, e -I /f(cr) values (i.e. v = —0.4) the time T for b to attain a 
given large value is longer than for v values corresponding to the larger x eh (<r) values 
(i.e. v = —1.4). In summary, the numerical calculations indicate that, for all 0 < v < 2, the 
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amplitude B(t ) will oscillate, with its magnitude rising sharply after a few oscillations; our 
results indicate that no ‘finite-time’ singularity of the evolution equation (5.6) occurs, instead 
the oscillations become increasingly wilder. 


5.3 Discussion 


We begin this discussion by continuing to consider the evolution equation (5.6) with non- 
linearity (solely) due to J\\ later in this subsection we consider the other (base) evolution 
equations. In the last section we reported that our numerical computations show that K(o) is 
everywhere negative (0 < a < 1) for all —2<v< 0. It was remarked that from these facts 
(and the magnitude of K = K(cr‘,u)) it is possible to deduce all the qualitative results of the 
solution properties of the evolution equation (5.6), without needing to perform the numerical 
integration. To illustrate that this is so, let us consider the right-hand side of equation (5.6). 
Initially b(T) > 0 and thus the right-hand side of equation is clearly negative (as K{a) < 0); 
thus this equation tells us that b(T ) will decrease in value until such time when the right-hand 
side of equation (5.6) is positive. The latter condition cannot be reached until b(T) becomes 
negative. However as the right-hand side of equation depends on all previous b(T ) values, there 
is a delay until the right-hand side of equation is actually positive; during this period b(T) has 
becomes more negative. Once the right-hand side of equation is positive, b(T) grows until it 
becomes positive. Again there is a delay until the right-hand side of equation is negative, at 
which time b(T) starts to decrease again, and so on. It is clear that the larger the typical 
magnitude of the overall kernel x 2 ~ 4t/ e~ x K(a), the more effect we can expect the nonlinearity 
to have i.e. larger oscillations at earlier T-values. Note that we would not expect a singularity 
to develop because of the smooth mechanism underlying the behavior of 6(T), as described 
above. In fact, it is possible to show analytically that the solution of the evolution equation 
(5.6) can only develop a singularity if /J K(a)da > 0 (see Churilov & Shukhman, 1988). 

Thus, disturbances initially governed by the evolution equation (5.6) would soon become 
so large in magnitude that their evolution would move into its second stage (region IIIc of 
Figure 5.1a) where the largest (leading order) term in the composite evolution equation (5.1) 
is the quintic nonlinearity. We note that this result is not at all surprising and could have been 
deduced as soon as the scales and terms of the critical-layer expansions (4.4a-d) were deduced: 
as viscosity does not enter our analysis at leading order, the base evolution equations due to 
Ji -nonlinearity, the OCF-nonlinearity and the quintic-nonlinearity will all lead to unbounded 
amplitude growth, whether by increasing disturbance oscillations or by a singularity occurring 
(they do not permit so-called equilibrium solutions). As viscosity does not enter at leading 
order it certainly cannot damp out the amplitude growth. We note that even with stronger 
viscosity effects, equilibrium states may still not necessarily be reached (see Goldstein &; Leib, 
1989). Moreover, we note that for the marginal instability case considered by Churilov & 
Shukhman (1988), the nonlinearity due to viscosity permitted singular solutions in which the 
amplitude became unbounded at a finite T. 
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6 Conclusion 


In this paper we have considered the nonlinear development of unstable disturbances in 
stratified shear flows where the Richardson number J is less than one quarter in value. Although 
such modes are initially fast growing, we have assumed that viscous-spreading effects result 
in them evolving in a linear fashion until they reach a state where their amplitudes are large 
enough but their growth rate have diminished significantly so that amplitude equations are 
derived using weakly nonlinear and unsteady critical layer theories. We have found that four 
different (base) integro-differential amplitude equations are possible, including one due to a 
novel mechanism (the l J\~ cubic nonlinearity). The relevant choice of amplitude equation, at 
a particular instance, being dependent on the relative sizes of the disturbance amplitude, the 
growth rate of the disturbance, its wavenumber and the viscosity of the fluid. This richness of 
choice of possible nonlinearities arises mathematically from the indicial Frobenius roots of the 
governing linear inviscid equation (the Taylor-Goldstein equation) not, in general, differing by 
an integer. The initial nonlinear evolution of a mode will be governed by an integro-differential 
amplitude equations with a cubic nonlinearity but the resulting significant increase in the size 
of the disturbance’s amplitude leads on to the next stage of the evolution process where the 
evolution of the mode is governed by an integro-differential amplitude equations with a quintic 
nonlinearity. Continued growth of the disturbance amplitude is expected during this stage, 
resulting in the effects of nonlinearity spreading to outside the critical level, by which time the 
flow has become fully nonlinear. 

We finish by mentioning some further points that may be worthy of investigation; mainly 
these are related to the relaxation of the assumptions made here in this paper. Obviously, 
the inclusion of viscosity at leading order in the critical layer (i.e. treating A h ~ 0(f)) would 
answer the questions posed in the previous section as to whether it can damp down the unlimited 
growth in the magnitude of the disturbances, allowing equilibrium solutions to exist. Another 
obvious extension of this work is to three-dimensional disturbances; then it may be possible 
to consider interactions between two or more disturbances. Comparision of our theory with 
experiments or large-scale numerical simulations is necessary at some point. As mentioned 
in the introduction, as well as geo-physical flow applications, the theory has applications to 
aerodynamical flow situations. In fact, in a related paper Blackaby, Dando & Hall (1993) 
apply the ideas contained in this paper to the problem concerning the nonlinear evolution of 
inviscid Gortler vortices in a three-dimensional boundary-layer; their problem is more complex 
than that for the model stratified flow considered here, as the longitudinal vortex problem 
is necessarily three-dimensional. The theory developed in this paper can also be applied to 
a nonlinear study of the inviscid vortex instabilities in the three-dimensional boundary-layer 
flow above a heated plate. 
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Appendix A 


This appendix contains some additional results, concerning the outer inviscid problem for 
the second harmonic, needed for the case v > 0. It is merely a generalisation of the unpublished 
Appendix B of Churilov & Shukhman (1987) to the case v / 0. 

Outside the critical layer, the first term of the second harmonic, satisfies 


L 2 ip 2 ] = Qi = Bl sinh 2l/+1 \y\ 


2 3Jo 

| I 

cosh 4 y 2 sinh 4 y j 


04 . 1 ) 


We introduce the functions f a and f b which satisfy the homogeneous equation L 2 rf = 0, such 
that 


f a = e -W*+3M(i + 0(e _2|y| )), f b = e * N/ * 5 + 5|! ' l (l + 0(e 2|vl )), as y -> ±oo. 


The solution of (A.l) is then written in the form 

4 1) = 5i(/(y) + C' ± / 0 (y)), (A.2) 

where C± are constants and / is a particular solution of equation (A.l) that does not contain 
e ±fc\A 2 +3|y| as y _► ioo. As y —* iO, 

/ = 2 (1 3 + 2 2 ‘' |/ ) ) |i>|- 1 ^ , '(l + 0{y 2 )) + + 0(y)) + + 0(y)y, (.4.3) 

where a and are to be determined. Also 

/.,6~a.,4!/l' + ‘'(l + 0(!/)) + A ll |!<l } - , '(l+0(!')). “ V - ±0. (A4) 
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where the <x a ,b stud j3 a ^ are fully determined by the definition of f a i.e. they are determined 
solely from the outer, inviscid problem. Thus, we can write 


= »r(l + 0(y 2 )) + Bl + A.±lw|*— ) (1 + 0 ( »)), (A.5a) 

o — J 


as y — ► ±0, where 

g 2 ± = a + C±oc a , h 2 ± = P + C±p a - 
It is easy to show from these relations that 

h 2+ - (6 2+ /6 2 _)(a/? 0 - P<**) aa[{g2+Ig2 ) _ (* 2+ /fc 2 _)]- 

To evaluate a and /?, we subtract the singular term in (A.3) by introducing 

/(y) = / - 2(3-1^ sinh_1+2,/ lyl cosh_2 y » 


(A. 56, c) 


(A.5d) 


(A. 6a) 


so that 

/ = a|j/|2 +,/ + /%|2 _I, + ---, as y-*± 0 , (A. 66) 

and f(y) satisfies 

L 2 f = Q 2 = 2(3 sinh_1+2,/ M cosh ” 2 y ( 8 tanh2 y ~( 1 + 2u )) ■ (A. 6c) 

Then 

/ = VMM*) - M*)Mv )] . ( A V 

where the Wronskian VK = f a fb — fafb — 2i/(ar a /5fc — a bPa) = —2 ky/k? + 3. Considering j/ — > 0, 
we see that 

a = Jo ~ aa ^ dz and P ~ Jo W^ b ^ a ~ ^ a ^ dz ' (^- 8 ) 


Appendix B 


Again following Churilov & Shukhman (1988), in accordance with (4.34), we write 

M+-M. = (M+ - M_) 11 + (M + - M — ) 02 + (M + - M_)si, (B.l) 

and after extensive calculations we obtain 
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(M+—M-)n = 


i\- 7v k$-* v 2-%+‘ / (l + 2v)T 2 {\ + t/) f°° , f 00 


* 2 r a - u) 


f dt f dt 4 f dt x f dt 2 

Jo Jo Jc Jc 


j dt 3 B(r -t-u- t x )B(r -t-t 4 - t 2 )B{T + t + t 1 + t 2 + t 4 - U)B{t -t-t 4 
J c 

{tit2t 3 )~*~' / t 4 + ‘'(h -t-h-ti- 2t 4 )-2-‘'(t 3 -t-t x - t 2 )->- u (t 3 -t x - t 2 )i 
(ts (Wi - h) 2 + t 2 (t 2 - t 3 ) 2 + ht 2 {2t 3 -t x - * 2 )] - i + h)-*-" x 

(^1 — ^2) 2_1/ (^l + ^2 + ^ 3 ) (2^3 — ^1 — ^2 )F — -V,~- 
H(t x - t 2 )H(t 3 -t-h-t 2- 2t 4 ), 


1 H.l 

2’ 4 2 't\)J 


wx 2i +v u~ 5u k^~ 3l/ (l -\-2t/)T(\ + v) f°° [°° f M 

(M+-M_) 02- ^(i -T) Jo d Jo dtl j c dt2 

J dt 3 B(r — t — t x )B(r — t — t 2 )B(r + t + t x + t 2 — t 3 )B(r — t — t 3 )(tt x )i + 
-2 t-t x - t 2 )~ 2 -'(t 1 - t 2 )^- v {t l + t 2 )- 2 - v { 2 t + t x + 1 2 ) x 
(< + *1 + (j)f(-j - - J - 1 ; J - §: |)«(‘. - - a - (1 - <2), 


(M + -M-) 31 = 


2-2 +l 't2" 5,/ A:2 _3,/ (l + 2i/) 2 r(| + I/) f°° f°° f°° 


*r 2 (i - u) 


tCC too too f 

/ dt dt\ / dt 3 / dt 2 
Jo Jo Jo Jc 


/ dt 4 B(r -t-t 3 - ti)B(r -t-t 3 - h )B(t + 2 t + t x + t 2 + t 3 - U)B{r ~ h 
J c 

{t\t 3 )^ v t 2 2 v t 4 5 ^(<4 — 3^ t x — t 2 — 2t 3 )~i~ v (ti + t 2 ) 2 l '(t x — t 2 y v 

(<4 - 3t - 2* 3 )* +l/ (* 4 - 2t 3 y?- l '(t 4 -St- 2 ti - 2 1 2 - 2t 3 )(2t 4 -6t- 3 t x - 3 1 2 - 

(t\ - <i)f(-5 - - j - 5; | - j: |)«(‘i - ‘»W‘4 -3 t-u-h- a,). 


- *a)x 
+^ x 


(5.2) 


X 

(B.3) 
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- 4f 3 ) x 
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Figure 5.1a A diagram of the various regimes of the critical layer for v < 0. ! 

I: viscous, steady critical layer; Landau-Stuart-Watson equation. 

II: strongly nonlinear, equilibrium critical layer; Benney & Bergeron theory. 

Ilia: unsteady critical layer; largest term in integro-differential equation (IDE) is cubic and 
due to viscosity. 

Illb: unsteady critical layer; largest term in IDE is cubic and due to J\. 

I lie: unsteady critical layer; largest term in IDE is quintic. 

The thick line on the diagram indicates the expected evolutionary path of the disturbance. : 
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Figure 5.1b A diagram of the various regimes of the critical layer for v > 0. 

I: viscous, steady critical layer; Landau-Stuart-Watson equation. 

II: strongly nonlinear, equilibrium critical layer; Benney & Bergeron theory. 

Ilia : unsteady critical layer; largest term in IDE is cubic and due to viscosity. 

Illb: unsteady critical layer; largest term in IDE is cubic and due to OCF. 

IIIc: unsteady critical layer; largest term in IDE is quintic. 

The thick line on the diagram indicates the expected evolutionary path of the disturbance. 
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Figure 5.4b Solution of the J\ evolution equation for v — —1.4. 
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